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Abstract. We study the stability of the results of 3-v oscillation analysis of atmospheric and reactor 
neutrino data under departures of the one-dominant mass scale approximation. In order to do so we 
perform the analysis of atmospheric and reactor neutrino data in terms of three-neutrino oscillations 
where the effect of both mass differences is explicitly considered. We study the allowed parameter space 
resulting from this analysis as a function of the mass splitting hierarchy parameter a = Am /AM 2 which 
parametrizes the departure from the one-dominant mass scale approximation. We consider schemes with 
both direct and inverted mass ordering. Our results show that in the analysis of atmospheric data the 
derived range of the largest mass splitting, AM 2 , is stable while the allowed ranges of mixing angles sin 2 #23 
and sin 2 #13 are wider than those obtained in the one-dominant mass scale approximation. Inclusion of 
the CHOOZ reactor data in the analysis results into the reduction of the parameter space in particular for 
the mixing angles. As a consequence the final allowed ranges of parameters from the combined analysis 
are only slightly broader than when obtained in the one-dominant mass scale approximation. 
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1 Introduction where cy = cos fly and sy = sin fly. Thus the pa- 

rameter set relevant for the joint study of solar and at- 
Super-Kamiokandc (SK) high statistics data § indicate mospheric conversions becomes six-dimensional: two mass 
that the observed deficit in the /i-like atmospheric events is differences, three mixing angles and one CP phase, 
due to the neutrinos arriving in the detector at large zenith Results from the analysis of solar and atmospheric data 



angles, strongly suggestive of the oscillation hypothe- in the framework of two-neutrino oscillation |15|,|16|,y,|]J 

sis. Similarly, the latest SNO results §§ in combination im P!y tnat thc required mass differences satisfy 
with the Super-Kamiokande data on the zenith angle de- Am 2 . <C Am 2 t (2) 

pendence and recoil energy spectrum of solar neutrinos [T 



and the Homestake f§, SAGE @, and GALLEX+GNO § 
^] experiments, have put on a firm observational basis the 
long-standing problem of solar neutrinos |L0] , strongly in- 
dicating the need for v e conversions. 



For sufficiently small Atuq the three-neutrino oscillation 
analysis of the atmospheric neutrino data can be per- 
formed in the one mass scale dominance approximation 
neglecting the effect of Amf^. In this approximation it 
follows that the atmospheric data analysis restricts three 



Altogether, thc solar and atmospheric neutrino anoma- r , , ,. , , „ o A o n 

° ' ^ of thc oscillation parameters, namely, Am^ — Am^ 2 , #23 

lies constitute the only solid present-day evidence for physics &nd ^ Thig ^ ^ approximation used in Ref . 

beyond the Standard Model [fill . It is clear that thc mini- ^ , c , , , . . , . , , a , r 

J 1 — ' Conversely for the solar neutrino analysis the effect of os- 

mum joint description of both anomalies requires neutrino ,. .,1 2 1. ^ 1 u j j 

■>t- 1 dilations with Am &tm can be taken to be averaged and 

conversions among all three known neutrinos. In the sim- , , , Air, j n rfmirm t ^ ■ 

° solar data constrains Am^vYi and fli3 ||19Lk1I. In this 



plest case of oscillations the latter are determined by the , . , , . , , £ r^r\r\r? 

approximation the reactor neutrino data from OHOUZ 

structure of the lepton mixing matrix [[L2l, which, in ad- . , • c , . , , ,. a 

1 Lt - 1 ' provides information on thc atmospheric mass dittcrcncc 

dition to thc Dirac-typc phase analogous to that of the , , , . . , n j j_i r-m 1 r u 

J ^ ^ b and the mixing angle fli3, and the CP phase becomes 

quark sector, contains two physical phases associated to , , , 

^ ' unobservablc. 

the Majorana character of neutrinos, which however are TT <• 1 j 

•' ' However the assumption of one mass scale dominance 

not relevant for neutrino oscillation 11131 and will be set to , , , . ., , r 

1 — v may not be a good approximation neither for reactor nor 

zero in what follows. In this case the mixing matrix U can r , , . , , . 1 r ^ 2 ■ 

for atmospheric data, m particular tor Zittiq in its up- 

be conveniently chosen in the form fll3| ln , , c ^ -\ r ^ 

J 1— t per allowed values. Effects of the departure of the one 

/ _ ■ \ mass scale dominance approximation in the analysis of the 

C13C12 S12C13 si 3 e 10 



-S12C23 - S23S13C12 e lS C23C12 - S23S13S12 e lS S23C13 
S23S12 - S13C23C12 e lS -S23C12 - S13S12C23 e lS C23C13 



CHOOZ reactor data 22 has been included in Ref. [ fL9| 
^||,^4|. For atmospheric neutrinos in Refs. [^,|fj,|2^] it 

(I) was shown that oscillations with two mass scales of the 



/ 
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order of could give a good description of the exist- the analysis including only the atmospheric data or the 
ing data for some specific values of the parameters. Some reactor data. 

analytical approximate expressions for the effects of keep- The outline of this paper is as follows. In Sec. [2] we 
ing both mass scales in the description of atmospheric describe our notation for the parameters relevant for at- 
neutrinos are presented in Rcfs. |^5|,^6,28 . Furthermore mospheric and reactor neutrino oscillations with two mass 
Refs. |^5|,^6| describe how the presence of the second mass scales and discuss the results for the relevant probabili- 
scale can lead to an increase in the number of sub-GeV ties. In Sec. || and ^ we show our results for the analysis 
electron events which seems to improve the description of of atmospheric neutrino and reactor data respectively. For 
the observed distribution. atmospheric neutrinos we include in our analysis all the 

contained events from the latest 1489 SK data set ||, as 

To further explore this possibility and to verify the 

well as the upward-going neutrino-induced muon fluxes 

consistency of the one-dominant mass scale approxima- 

from both SK and the MACRO detector The results 

tion we present in this work the result of the analysis of 

for the combined analysis are described in Sec. H Finally 

the atmospheric and reactor neutrino data in terms of 

m Sec. | we summarize the work and present our conclu- 

three-neutrino oscillations where the effect of both mass 

sions. 

differences is explicitly considered and we compare our 
results with those obtained under the assumption of one- 

, . , , ^ . . , , , , , , , 2 Three Neutrino Oscillations with Two Mass 

dominant scale. Our aim is to study now/wncthcr the al- 
lowed parameter space is modified as a function of the Scales 
ratio between the two mass scales. Our study allow us to 

In this section we review the theoretical calculation of the 

establish the stability of the derived ranges of parameters 

conversion probabilities for atmospheric and reactor neu- 

for the large mass scale and mixings #23 and #13 inde- 

trinos in the framework of three-neutrino mixing, in order 

pendently of the exact value of the solar small scale and 

to set our notation and to clarify the approximations used 

mixing #12 for which we only chose it to be within the 

in the evaluation of such probabilities. 

favoured LMA region. Our results show that the allowed 

In general, the determination of the oscillation prob- 

ranges of parameters from the combined atmospheric plus 

abilities for atmospheric neutrinos require the solution of 

reactor data analysis are only slightly broader than when 

the Schrodinger evolution equation of the neutrino sys- 

obtained in the one-dominant mass scale approximation. 

tern in the Earth-matter background. For a three-flavour 

Thus our main conclusion is that the approximation is self- 
scenario, this equation reads 

consistent. To establish the relevance of each data sample 

on this conclusion we also present the partial results of = Hz/, H = U • Hq ■ + V , (3) 
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; this case we can have the two mass ordering: 



NORMAL 



Normal: AM 2 = Am\ x = m\ 



'21 



Inverted: AM 2 = —Am 



INVERTED 



Am 2 = Am\ x 



2 2 
77^2 — m\ . 



(6) 



(7) 



We define the mass splitting hierarchy parameter 



Fig. 1. Our convention for the mass splitting and ordering. 



(8) 



where U is the unitary matrix connecting the flavour basis AM 2 
and the mass basis in vacuum and which can be parametrized which parametrizes the departure from the one-dominant 
as in Eq. (jj]). On the other hand H$ and V are given as mass sca le approximation in the analysis of atmospheric 



and reactor neutrinos. 

In this convention, for both Normal or Inverted schemes, 
the mixing angles in Eq. (Q) are such that in the one mass 
where v = {y ei v^v T ). We have denoted by H„ the vac- dominance approximation in which AM 2 {Am 2 ) deter- 



H ° = dr diag M m !i> 4m 3i) 



V = diag (±V2G F N e , 0, 



(4) 
(5) 



uum Hamiltonian, while V describes charged-current for- mines the oscillation length of atmospheric (solar) ncu- 

ward interactions in matter In Eq. (§), the sign + trinos, 9 2 3 is the mixing angle relevant for atmospheric 

(-) refers to neutrinos (antineutrinos) , G F is the Fermi oscillations while 6»i 2 is the relevant one for solar oscilla- 

coupling constant and N e is electron number density in tions, and 9i 3 is mostly constrained by reactor data. In 

the Earth (note also that for antineutrinos, the phase S the likely situation in which the solar solution is LMA, 



has to be replaced with —5). 

The angles 9ij can be taken without any loss of gen- 



?i2 is mainly restricted to lie in the first octant. 

We will restrict ourselves to the CP conserving sce- 



crality to lie in the first quadrant 0y 6 [0, 7r/2]. Con- nario. CP conservation implies that the lepton phase S is 
cerning the CP violating phase 5 we chose the convention either zero or ir Q . As we will see, for non- vanishing a 
< S < 7T and two choices of mass ordering (See Fig. 0) and 6*i 3 the analysis of atmospheric neutrinos is not ex- 
one with mi < rri2 < m 3 which we will denote as Normal actly the same for these two possible CP conserving values 
and other with m 3 < mi < m 2 which we will denote as of S and we characterize these two possibilities in terms of 
Inverted (for a recent discussion on other conventions see, cos<5 = ±1. 

for instance [p2|). We define as AM 2 > the large mass For a = 813 = 0, atmospheric neutrinos involve only 



splitting in the problem and Am 2 > the small one. In Vn — * v T conversions, and in this case there are no matter 
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effects, so that the solution of Eq. (||) is straightforward 
and the conversion probability takes the well-known vac- 
uum form 



= 1 - sin 2 (20 23 ) sin 1 



/ AM 2 L 
{ 4E» 



(9) 



where L is the path-length traveled by neutrinos of energy 
E v . 

On the other hand, in the general case of three- neutrino 
scenario with #13 ^ or a 7^ the presence of the mat- 
ter potentials become relevant and it requires a numer- 
ical solution of the evolution equations in order to ob- 
tain the oscillation probabilities for atmospheric neutrinos 
P a f3, which are different for neutrinos and anti-neutrinos 
because of the reversal of sign in Eq. (|^). In our calcula- 
tions, we use for the matter density profile of the Earth the 
approximate analytic parametrization given in Ref. J34| of 
the PREM model of the Earth 




Standard Model 
a=0.00, s, 3 =0.0, S23=0.35 
a=0.00, sf 3=0. 1 , S2 3 =0.35 
a=0.00, sf 3 =0.1, S2 3 =0.65 
a=0.15, s^ 3 =0.0, S23=0.35 
a=0.15, s?,=0.0, sL=0.65 



-1 -0.8 -0.6 -0.4 -0.2 
COS 6 



In Figs. ^ and || we plot the angular distribution of 
atmospheric v e and for non-vanishing values of a or 

013 obtained from our numerical calculations. As seen in Fig _ 2 Zenith _ angle distr ibutions (normalized to the no- 

these figures the main effect of a small but non- vanishing osci u a tion prediction) for the Super-Kamiokande e-like and 

a is mostly observable for sub-GeV electrons, although nke contained events, for the Super-Kamiokande stopping and 

some effect is also visible for multi-GeV electrons and sub- through-going muon events and for Macro upgoing muons. The 

GeV muons, and it can result either in an increase or in a various dashed lines are the expected distributions for the Nor- 

decrease of the expected number of events with respect to mal mass ordering with AM 2 = 3 x 10" 3 eV 2 , tan 2 8 12 = 0.45 

the a = prediction depending on whether 6» 23 is in the and several values of sin 2 X3 and sin 2 6» 23 as given in the figure, 
first or second octant. This behaviour can be understood 
in terms of the approximate analytical expressions. For 
instance for #13 = we find (in agreement with the results 
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0.8 



SK multi-GeV (e) 

-0.5 0.5 



1 a4 -1 

0. 




SK multi-GeV (u) 

-0.5 0.5 



SK stop (u) 

1 -0.8 -0.6 -0.4 -0.2 
COS 6 

Inverted 



Standard Model 
a=0.00, s* 3 =0.0, s|j=0.35 
a=0.00, S, 3 =0.1, s|j=0.35 
a=0.00, S%=OA, s|j=0.65 
a=0.15, s^ 3 =0.0, S23=0.35 
a=0.15, S?,=0.0, sL=0.65 



-1 -0.8 -0.6 -0.4 -0.2 
COS 9 



Fig. 3. Same as Fig. g but for Inverted mass ordering, 
in Ref. §|) 



N e0 



1 = P e2 f{c: 



1 



N„-Nja = 0) 



N, 



/,0 



23 



:) 



(10) 

(11) 



the expected number of muon-like events for a = and 
P e 2 is the dominant a-dependent term in the probabilities, 
averaged over energy and zenith angle. For neutrinos we 
have: 



sin 2 2f5i2 m sin 2 



sin26»i2, m 



Am 2 L sin 26»i2 

AE V sin 26»i2, m 
sin 26>i2 



(12) 



y/{cos 26*12 - 2E v V e /Am 2 ) 2 + (sin26ii2 



which for Am 2 <C 1E v V e reduces to: 



P e2 = a 2 sin 2 26*12 



(AM 2 
\2EV e 



sin — . (13) 



factor c 2 3 



<-2.i 



According to Eqs. ( p^| ) and ( pi] ) the sign of the shift in the 
number of predicted events with respect to the results in 
the one mass scale dominance approximation is opposite 
for electron and muon-like events and it depends on the 
0.5. So for 6*23 in the first octant, 
c 23 > O-^) there is an increase (decrease) in the number 
of electron (muon) events as compared to the a = case. 
For #23 in the second octant the opposite holds. We also 
see that the net shift is larger for electron events than 
for muon events by a factor c 2 3 /f. Notice that, despite 
Eq. ( |T3| ) looks order a 2 , its numerical value for sub-GeV 
electrons is large due to the factor AM 2 / (2EV e ) as can 
be seen from the figures. At higher energies, for up-going 
muons the effect is negligible. 

For the sake of comparison we also show in the figures 
the behaviour with non-vanishing value of 6*13 in the one 



where N e o and N^q are the expected number of electron mass scale dominance approximation. As seen in the figure 

and muon-like events in the absence of oscillations in the the effect is most important for the electron events and can 

relevant energy and angular bin and f = N^/NeQ. For be understood as follows. For the case of constant matter 

instance, for sub-GeV events f ~ 2. Here JV M (a = 0) is density the expected flux of v e events in the one mass scale 
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dominance approximation we find whether the Am 2 oscillations are above (Am 2 cos 2(9i2 > 

i 2E u V e ) of below (Am 2 cos 28 \i < 2E l/ V e ) the resonance. 
^_l = P eM f(4- F ) (14) 

e0 For very small a (Am 2 <C 2E„V e ) the interference term 

where is proportional to a and it also changes sign for neutrinos 

/ AM 2 L sin 26*i •? \ and antineutrinos. In summary the effect of non- vanishing 

Pen = 4s? 3 ,mCl3 >m ^ ( , P „ ; „ ^ ) > (15) 



sin20i 3 . m 



V 4E„ sin2(9i3„ 
sin 2^i3 



?i3 and a in the expected number and distribution of at- 



^(cos 26*13 T 2E y V e / AM 2 ) 2 + (sin26>i3) 2 mospheric neutrino events can have opposite signs and this 

and the - (+) sign applies for the Normal (Inverted) case can lead to a P artial cancellation between both contribu- 

(similar expression is presented, for instance, in the last tions - This results into a loss of sensitivity of the analysis 

article in Ref. [||| and in Ref. So for 6» 23 in the first to both parameters. 

octant (s| 3 < 0.5) there is a decrease in the number of To a nalyze the CHOOZ constraints we need to evalu- 

electron events as compared to the 13 case. For sub-GeV ate the survival probability for D e of average energy E ~ 



events, the matter term in Eq. © can be neglected and few MeV at a distance of L ~ 1 Km. For these values of 

the effect of a non-vanishing 8 13 is the same for Normal energy and distance, one can safely neglect Earth matter 

and Inverted ordering. For multi-GeV and upgoing muon effects. The survival probability takes the analytical form: 

events matter effects start playing a role and the effect be- d chooz i „„„4/i • 2 /j 2 f Am 21 L 

r J r = 1 — COS f 13 sin It) 12 sin — — — 

comes slightly larger for the Normal case where the matter 



sin 2 26» 



13 



enhancement is in the neutrino channel. 



cos 2 On sm 2 ( ^£ 



• sin"* 6>i2 sur ' 



The situation becomes more involved when both a and " y AE„ 

( AM 2 L 



(17) 



#13 are different from zero. For instance, in lowest order in ~ 1 — sin 26>i3 sin 

a S13 the expected number of sub-GeV v e events is (after 
averaging the AM 2 L/E oscillations) 



V 4£„ 

where the second equality holds under the approxima- 
tion Am 2 <C Ey/L which can only be safely made for 
Am 2 < 3 x 1(T 4 cV 2 . Eq. @ is valid for both Normal 
and Inverted ordering with the identifications in Eq.(|^) 
^ m 2 Yj — s i n 28 \~ an< ^ ^' C '' ® respectively It results that the probability for 



~ 1 = Pe2 f (4 -h + P e „ r (s 2 23 - i) (16) 

iv e o T T 

+ — cos S sin 26* 13 sin 26* 2 3 sin 26 ) i2, m cos 26*12, m 

. 2 
sin 



4E„ sin 26*12,™ / Normal and Inverted schemes is the same with the ex- 
(this expression is in agreement with the results in Ref. |p6f ). change sin 2 6*12 <-> cos 2 6*12. Thus in general the analysis 
From this equation we see that the interference term (the of the CHOOZ reactor data involves four oscillation pa- 
third term in the right hand side) can have either sign rameters: AM 2 , 6*13, Am 2 , and 6*12- From Eq. ( |l7| ) we see 
depending on cos<5. It also changes sign depending on that for a given value of 6*12 and AM 2 the effect of a non- 
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vanishing value of either #13 or Am 2 is the decrease of the global analysis of solar neutrino data by choosing a char- 



survival probability. acteristic value tan 2 #i 2 = 0.45 [|J|,|T6|. We have repeated 

our analysis for different values of 812 and we have found 

3 Atmospheric Neutrino Analysis that the maximum effect due to the variation of 9i2 is a 

shift on A\ 2 ~ 1 and it is therefore unobservable. Further- 
In our statistical analysis of the atmospheric neutrino events 

more we have verified that the atmospheric data analysis 

we use all the samples of SK data: e-like and /j,-like sam- 

does not provide enough precision to test the possibility 

pies of sub- and multi-GeV § data, each given as a 10-bin 

of non- vanishing CP violation. 

zenith-angle distribution , and upgoing muon data includ- 

We first present the results of the allowed parameters 

ing the stopping (5 bins in zenith angle) and through-going 

for the global combination of atmospheric observables. 

(10 angular bins) muon fluxes. We have also included the 

__ Notice that since the parameter space we study is four- 

latest MACRO [BOf upgoing muon samples, with 10 angu- 

dimensional the allowed regions for a given CL are de- 

lar bins. So we have a total of 65 independent inputs. 

fined as the set of points satisfying the condition for four 

For details on the statistical analysis applied to the dif- 

degrees of freedom (d.o.f.) 

ferent observables, we refer to the first reference in Refs. |17| 

and [091 . As discussed in the previous section, the analysis 2 mi/2/i n a 2\ 2 ^ a a j c\ 

min — \ { ^1 J ] 4 U.O.I.J 

of the atmospheric neutrino data for three neutrino oscil- n oN 

[16) 

lations with two mass scales involves six parameters: two where 4 ±q f } = ? ^ g ^ ^ and ^ for 

mass differences, three mixing angles and one CP phase. CL = QQ% ^ ^ gg% and ggj3% _ ^ respectively) and 

Our aim is to study the modification on the resulting al- 9 . . n n , , . . , , r ,. . , 

Xatm min 1S ^ e global minimum 111 trie tour-dimensional 

lowed ranges of the parameters AM 2 , sin 2 $93 and sin 2 9\-i , , n , . , , , , „ ,, n 

01 ' space, lhe best fit point used to define the allowed pa- 
due to the deviations from the one-dominant mass scale . r n , , 

rameter space is found to be: 

approximation, i.e. for Am 2 ^ (or equivalently for non- 
vanishing values of mass splitting hierarchy parameter a). Z\M — 3.3 x 10 eV 
In what follows, for the sake of simplicity, we will restrict sin 2 6*23 = 0.46 

ourselves to the CP conserving scenario but we will distin- sm 2 q 13 — q (19) 

guish the two possible CP conserving values of 8 and we . o 1 n m-3 i n on\ 

v Am = 1.0 x 10 ev (a = 0.30) 

characterize these two possibilities in terms of cos 8 = ±1. 

v . = 3Q n 

Aatm,mm utj.\j 

We will show the results for Normal and Inverted schemes. 

Furthermore in most of our study we will keep the mix- (for 61 d.o.f.) and it corresponds to Normal ordering al- 
ing angle #12 to be within the LMA range favoured in the though the difference with the Inverted ordering (A\ = 
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(for 62 d.o.f., one more than in the a-unconstrained case) 
which is independent of the choice cos S = ±1 and corre- 
sponds to Inverted schemes. Notice that this is the min- 
imum used to obtain the allowed parameter space in the 
one-dominant mass scale approximation [the upper panels 
in Fig. |t] and the ranges in Eq. (p2[)]. since in this case we 
are fixing a priori a = 0. 

In summary, we find that allowing for a non-zero value 
of a very mildly improves the quality of the global fit. 
This result is driven by the better description of the sub- 
GeV data which is attainable for a non-zero a value, and 




0.2 0.3 
a = Am 2 / AM 2 

Fig. 4. Dependence of the A\lt m and ^Xatm+CHOOz functions 
on the mass splitting parameter a for tan 2 612 = 0.45, for the 
analysis of atmospheric neutrinos alone (lower panel) and also 
in combination with the CHOOZ reactor data (upper panel). 



drives the best fit point to the first octant of the mix- 
0.1) is not statistically significant []. The point given in in S an § le ^3 for which the expected number of sub-GeV 
Eq. © is the global minimum used in the construction electrons is larger as compared to the pure i/„ -» v T sce- 
of the A X l tm function shown in Figs. | and g of the al- nario ' as illustrated in Figs. | and |. We find, however, 
lowed parameter space shown in Fig. § and in the lower that the analysis of the atmospheric neutrino data does 



panels of Fig. [?], and of the ranges in Eq. (f2l|). 

This result can be compared with the best fit point 
obtained in the one-dominant mass scale approximation 
a = 

AM 2 = 3.0 x eV 2 
sin 2 6>23 = 0.54 



sin 013 



0.14 
39.6 



(20) 



1 The careful reader may notice that the \ 2 per d.o.f. seems 



not show a strong dependence on large a values. In the 
lower panel of Fig. ^ we show the dependence of Z\x 2 tm 
on a. In this plot all the neutrino oscillation parameters 
which are not displayed have been "integrated out" , i. e. 
the Axit m function is minimized with respect to all the 
non-displayed variables. From this figure we see that the 
fit to atmospheric neutrinos is only weakly sensitive to the 
value of a. 

In Fig. |^ we present sections of the allowed volume 
in the plane (cos S sin 2 623, AM 2 ) for different values of 



too good. This was already the case for the previous Su- sin 2 6 13 and for values of Am 2 = (first row) and Am 2 = 
perKamiokande data sample and it is partly due to the very 3 x 10~ 4 eV 2 (second row) which is the (maximum al- 
good agreement of the multi-GeV electron distributions with lowed value by the present analysis of solar neutrino data 



their no-oscillation expectations 



including the latest 1496 days of SK and and day-night 



10 
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Fig. 5. 90%, 95%, 99% and 3a (4 d.o.f.) allowed regions in the (sin 2 823, AM 2 ) plane, for different values of sin 2 13 and Am 2 , 
from the analysis of the atmospheric neutrino data. The global minimum used to define the allowed regions is given in Eq. (^); 
the local minima are marked with a dot. 
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spectrum of SNO data . For illustration we also show • For a — (no dependence on cos S) 
the corresponding regions for the "democratic" scenario 
a = 0.5 (Am 2 = AM 2 /2). We display the corresponding 
sections for the Normal and Inverted schemes. 



Normal Inverted 



1.3 < j^v, < 8.1 1.3 < j^Ly, < 10.0 



0.32 < sin 2 6*23 < 0.79 0.32 < sin 2 d 23 < 0.78 

Comparing the sections in Fig. || for a = with the sin 2 8 13 < 0.49 sin 2 6 13 < 0.58 

(22) 

corresponding sections for non vanishing a values we find 

that substantial differences appear although mainly for Comparing the ranges in Eqs. @ and @ we see 
large values of $13. However from these figures one also that the parameter which is less sensitive to the depar- 
realizes that even for large values of a the allowed region ture from the one mass scale dominance approximation 
does not extend to a very different range of AM 2 . Con- is AM 2 , while sin 2 6*13 is the mostly affected, in particu- 
versely, the mixing angles 623 and #13 can become less lar for the Inverted scheme for which no upper bound on 
constrained when the case a 7^ is considered. sin 2 #13 is derived from the analysis. The careful reader 

may notice that for the Normal ordering, the bound on 



To further quantify these effects we plot in Fig. |fj| the 
dependence of Ax 2 tm on AM 2 , 823 and #13, respectively, 
for different values of a, after minimizing with respect to 
all the non-displayed variables. From these figures we can 
read the 3er allowed ranges for the different parameters (1 
d.o.f.): 

• For arbitrary a 

Normal Inverted 
1-3 < W A ?% < 8.1 1-2 < W ^ M * V , < 9.6 

cos (5 = 1 

0.22 < sin 2 6*23 < 0.79 0.14 < sin 2 8 23 < 0.78 

sin 2 8 13 < 0.48 sin 2 6 13 < 0.58 

cos S = — 1 

0.19 < sin 2 6*23 < 0.79 0.22 < sin 2 9 23 < 0.95 

sin 2 6»i3 < 0.48 sin 2 13 < 1 

(21) 



613 for arbitrary a can be stronger than for a = 0. This 
is due to the fact that the ranges in Eqs. (|2|) and (|22 
are defined in terms of 3er shifts in the x 2 function with 
respect to the minima in Eqs. (|lj]) and (^0|) respectively 
[see explanation below Eq. (po])]. 

Finally we show in Fig. [?] the 2-dimensional allowed 
regions in (cos J sin 2 #23, AM 2 ) from the analysis of the 
atmospheric neutrino data independently of the values of 
a and #13. In constructing these regions for each value of 
AM 2 and cos S sin 2 #23 we have minimized on the oscilla- 
tion parameters Am 2 and #13 so the they are defined in 
terms of A X 2 for 2 d.o.f. (A X 2 = 4.61, 5.99, 9.21, 11.8 
for 90%, 95%, 99% CL and 3cr respectively). For the sake 
of comparison we show in the figure the corresponding re- 
gions for a = 0. From the figure we see that the differences 
are larger for the Inverted scheme. 
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Fig. 6. Dependence of the Axitm function on the mixing angles cos S sin 2 623 and cos S sin 2 #13 and on the large mass scale 
AM 2 , for different values of a and for the Normal (upper panels) and Inverted (lower panels) cases. See text for details. 



4 Analysis of CHOOZ Data 



The CHOOZ experiment [£2| searched for disappearance 
of v e produced in a power station with two pressurized- 
water nuclear reactors with a total thermal power of 8.5 GW 
(thermal). At the detector, located at L ~ 1 Km from the 
reactors, the v e reaction signature is the delayed coinci- 
dence between the prompt e + signal and the signal due to 
the neutron capture in the Gd- loaded scintillator. Their 
measured vs. expected ratio, averaged over the neutrino 
energy spectrum is 



R= 1.01 ± 2.8 %(stat) ±2.7 %(syst) 



(23) 



Thus no evidence was found for a deficit of measured vs. 
expected neutrino interactions, and they derive from the 
data exclusion plots in the plane of the oscillation param- 
eters (Am 2 , sin 2 29) in the simple two- neutrino oscillation 
scheme. At 90% CL they exclude the region given approx- 
imately by Am 2 > 7 ■ 10 eV 2 for maximum mixing, and 
by sin 2 (29) > 0.10 for large Am 2 . Similar searches have 
been performed at the Palo Verde Reactor Experiment 
[ p7[ leading to slightly weaker bounds. 

In order to combine the CHOOZ bound with the re- 
sults from our analysis of atmospheric neutrino data in 
the framework of three-neutrino mixing we have first per- 
formed our own analysis of the CHOOZ data. Using as ex- 
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Fig. 7. 90%, 95%, 99% and 3cr (2 d.o.f.) allowed regions in 
the (cos r5 sin 2 #23, AM 2 ) plane from the analysis of the atmo- 
spheric neutrino data, for the Normal (left panels) and Inverted 
(right panels) cases, for tan 2 612 = 0.45 and for arbitrary val- 



in Sec. g for the analysis of the reactor data the relevant 
oscillation probability depends in general on the four pa- 
rameters 6*12 , Am 2 , 613, and AM 2 . In Fig. || we show the 
excluded regions at 90, 95 and 99% CL and 3er in the 
(sin 2 # 13 , AM 2 ) plane from our analysis of the CHOOZ 
data for several values of Am 2 and tan 2 #12 = 0.45; for 
the sake of comparison with the 2-family analysis we have 
defined the allowed regions for 2 d.o.f. (^Xchooz = 4.61, 
5.99, 9.21, 11.83 respectively). In the left (right) panel we 
show the results for the Normal (Inverted) scheme. We see 
that the presence of a non-vanishing value of Am 2 results 
into a slightly smaller allowed range of {AM 2 , sin 2 0i 3 ). 
For the chosen value of tan 2 f5i 2 the reduction for smaller 
values of AM 2 is slightly more significant for the Normal 
than for the Inverted scheme as also shown in Ref. |E3|. 



ues of #13 and a (lower panels). See text for details. The upper 
panels correspond to the case a = 0. The best fit point in each This can be easil y understood from the expression of the 
case is marked with a star. The local minima are marked with survival probability: from Eq. @, we get 
a dot. 



perimental input their measured ratio ( p3[ ) p2[ and com- 
paring it with the theoretical expectations we define the 
Xchooz function. We verified that with our Xq H ooz func- 
tion and using the statistical criteria for two degrees of 
freedom we reproduce the excluded regions given in Ref. [^2| 
as can be seen in the upper row of Fig. ||0 As discussed 



pCHOOZ pCHOOZ _ • 2 ofl /„-.„2 n • 2 a \ 

-^ee.NOR ~ ^ee.INV -"Sin ZU 13 {COS U 12 - Sin V U ) 

( AM 2 L 



sin 



— sm 



{AM 2 - Am 2 )L 



(24) 



2 For the sake of simplicity we chose not to include the en- 
ergy dependence of the CHOOZ data which adds very little 



Thus for #i2 < \ the survival probability is smaller for the 
Normal ordering than for the Inverted one, which leads to 
the stronger constraint. For AMI 2 3> Am 2 , Eq. (|24|) van- 
ishes and the excluded regions in the two schemes become 
indistinguishable. 

5 Combined Analysis 



to the knowledge of the parameter space as can be seen by 

comparing our results in Fig. ||with those of the CHOOZ col- We now describe the effect of including the CHOOZ re- 

laboration Q or the corresponding ones in the analysis of actor data together with the atmospheric data samples 

Ref. in a combined 3- neutrino \ 2 analysis. The results of this 
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Fig. 8. 90%, 95%, 99% and 3cr (2 d.o.f.) allowed regions from 
the analysis of the CHOOZ reactor data in the (sin 2 #13, AM 2 ) 
plane for different values of Am 2 (tan 2 #12 = 0.45), for the 
Normal (left panels) and Inverted (right panels) cases. 



the point given in Eq. (|2^) is the global minimum used in 
the construction of the ^Xatm+CHOOZ function shown in 
Figs. [I| and ^|, of the allowed regions in the lower panels of 
Fig. [h], and in the ranges in Eq. (|27|). 
For a = the best fit point is at: 



AM 2 = 2.5 x 1(T 3 eV 2 
sin 2 6> 23 = 0.49 (~ 0.51) 
0.005 



sin 0i3 



(26) 



Xatm+CHOOZ.min — ^0.2 



for 63 d.o.f. This is the minimum used to obtained the 
allowed parameters in the one-dominant mass scale ap- 
proximation: the upper panels in Fig. |l^ and the ranges 
analysis are summarized in the upper panel of Fig. ^ and m jrjq (ps|), 

in Figs. As in Sec. |, in most of results shown here In the upper panel of Fig | we show the dependence of 

we fix the mixing angle tan 2 9 12 = 0.45 and study how A\; a tm+CHOOZ 011 a - From this fi S ure we see that thc in " 
the allowed ranges of the parameters AM 2 , sin 2 9 23 and c i usion f the CHOOZ reactor data results into a stronger 



sin 9i3 depend on a. 



dependence of the analysis on the value of Am 2 (or equiv- 



We find that, in this case, the best fit point for the a i en tly on a) and large values of thc mass splitting hierar- 
combined analysis of atmospheric and CHOOZ data is chy parame ter become disfavoured. Also the dependence 



practically insensitive to the choice of Normal or Inverted 
schemes and: 

AM 2 = 2.8 x 10~ 3 eV 2 
sin 2 6» 23 = 0.46 
sin 2 6>i 3 = 0. 

Am 2 = 2.8 x 10" 4 eV 2 (a = 0.1) 



(25) 



is stronger for the Normal scheme, as expected (see dis- 
cussion below Eq. (p^)). As a consequence the ranges of 
mixing parameters - which, in the analysis of atmospheric 
data alone, were broaden in presence of large values of a - 
are expected to become narrower with the inclusion of the 
CHOOZ data in the analysis. 

This effect is explicitly shown in Fig. H, where we plot 



Xatm+CHOOZ,: 



39.8 



the dependence of the ^Xatm+CHOOZ 011 AM 2 



and 



for 62 d.o.f. and cos<5 = ±1. Notice than in our analysis #13, respectively, for different values of a (to be compared 
the CHOOZ data adds only one data point. Note that with the corresponding Fig. || for the analysis of the at- 
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Fig. 9. Dependence of the ^Xatm+CHOOz function on the mixing angles cos S sin 2 #23 and cos S sin 2 #13 and on the large mass 
scale AM 2 , for different values of a and for the Inverted case. See text for details. 



mospheric data). Fig. [9] is shown for the Inverted scheme, 
(the corresponding figure for the Normal scheme is very 
similar). The figure illustrates that indeed the inclusion of 
the CHOOZ data in the analysis results into a reduction 
of the allowed ranges for the mixing angles, in particular 
613. From this analysis we obtain the 3ct allowed (1 d.o.f.) 
bounds: 

• For arbitrary a 



Normal 



1.3 < 



AMr 

10 -3 c y2 



< 5.4 



Inverted 



1.0 v 10 _ 3 cV 2 v o.z 



cos 5 = ±1 



0.26 < sin 2 6*23 < 0.71 
sin 2 6>i 3 < 0.06 

• For a = 



Normal 
1 ^ < am2 < c; 1 

0.31 < sin 2 d 23 < 0.71 
sin 2 6>i 3 < 0.07 



0.26 < sin 2 023 < 0.70 



sin 2 0i3 < 0.07 



(27) 



Inverted 
1.3 < < 5.0 

0.31 < sin 2 23 < 0.70 
sin 2 0i3 < 0.07 



(28) 



Comparing with the results in Eq. (|2lj) we see that includ- 
ing the CHOOZ reactor data reduces the effect on the final 
allowed range of parameters arising from allowing depar- 
tures from the one mass scale dominance approximation. 
In other words the ranges in Eq. (|27]) and ( p8| ) are not 
very different. 

Fig. [To] shows the global 2-dimensional allowed regions 
in (cos S sin 2 023, A M 2 ) from the analysis of the atmo- 
spheric neutrino and CHOOZ reactor data for optimized 
values a and 0i3 as well as the results for the one mass 
scale dominance approximation a = case. Comparison 
with Fig. ^ shows that after including the CHOOZ reactor 
data the allowed range of parameters AM 2 and sin 2 023 
becomes more "robust" and it is almost independent of 
the Normal or Inverted ordering of the masses or of the 
particular choice of cos (5 = ±1. 

How large would have a and/or 0i3 to be for Z-v effects 
to be visible in the global analysis?. We find that in order 
to have a 3tr effect on the global analysis either tan 2 0i3 
should be larger than 0.07 or a should be larger than 0.4. 
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Fig. 10. 90%, 95%, 99% and 3a (2 d.o.f.) allowed regions in 
the (cos 5 sin 2 823, AM 2 ) plane, from the analysis of the atmo- 
spheric and CHOOZ neutrino data with tan 2 612 = 0.45, for 
the Normal (left panels) and Inverted (right panels) cases and mospheric data alone or in combination with the CHOOZ 



have studied the effect of keeping the two mass scales on 
the combined analysis. 

In general the analysis of atmospheric data involves 
six parameters: two mass differences, which we denote by 
AM 2 and Am 2 , three mixing angles (823, 9 13 and f5 12 ) 
and one CP phase (5). The analysis of the reactor data 
involves four of these parameters, namely, AM 2 and 5m 2 , 
#13 and 0\2- For the sake of simplicity we have concen- 
trated on the dependence on Am 2 while keeping the CP 
phase fixed to CP conserving values and the mixing angle 
#12 to be within the LMA range favoured in the global 
analysis of solar neutrino data by choosing a character- 
istic value tan 2 6*12 = 0.45. We have verified that the at- 



for arbitrary values of 6> 13 and a (lower panels). See text for data is not sensitive enough to give any constraint on the 
details. The upper panels correspond to the case a = 0. The possibility of CP violation nor to variations of the tan 2 6 12 
best fit point in each case is marked with a star. The local within the allowed LMA range. Thus our conclusions are 



minima are marked with a dot. 



6 Summary 



In this article we have explored the effect of keeping the 
two mass scales on the three-flavour oscillation analysis of 
the atmospheric and reactor neutrino data. First we have 
performed the independent analyses of the atmospheric 
neutrino data and of the CHOOZ data. We have studied 
the allowed parameter space resulting from these analyses 
as a function of the mass splitting hierarchy parameter 



robust. 

Our results can be summarized as follows: 

— The dominant effect of a non- vanishing value of a in 
the atmospheric neutrino events is an increase (de- 
crease) of the expected number of contained 24 for 623 
in the first (second octant) as previously discussed in 

Refs. mm. 

— In the predicted atmospheric neutrino events the ef- 
fects of a non- vanishing a and of the mixing angle #13 
can have opposite signs and certain degree of cancel- 
lation may occur between both effects. 



a = Am 2 1 AM 2 which parametrizes the departure from — The survival probability of v e at CHOOZ decreases 



the one-dominant mass scale approximation. Finally we 



for increasing values of #13 and a, so that the effect 
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of both parameters is additive in the CHOOZ reactor 
data. For 6*12 < f the effect of Am 2 is slightly stronger 



for the Normal mass ordering p3,24 



— Allowing for a non-zero value of a very mildly improves 
the quality of the atmospheric neutrino fit as a conse- 
quence of the better description of the sub-GeV elec- 
tron data. This effect drives the best fit point to the 
first octant of the mixing angle 623- 

— Still the fit to atmospheric neutrinos is very insensitive 
to large values of a as long as all other parameters are 
allowed to vary accordingly. 

— As a consequence the allowed range of sin 2 6*13 and 
sin 2 023 from the atmospheric neutrino data analysis 
becomes, in general, broader than the one for the a — 
case. 

— On the other hand the allowed range of AM 2 obtained 
from the atmospheric neutrino data fit is stable under 
departures from the one mass scale dominance approx- 
imation. 

— The inclusion of the CHOOZ reactor data in the analy- 
sis leads to a stronger dependence of the results on the 
value of a, with smaller values of a and #13 favoured. 

— As a consequence the final determination of the al- 
lowed ranges for both AM 2 and the mixing angles #23 
and 0i3 is very robust and the ranges are only slightly 
different from those obtained in the one mass scale 
dominance approximation. 
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